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The Infrared Array Camera (IRAC) on the Spitzer Space Telescope is abso- 
lutely calibrated by comparing photometry on a set of A stars near the north 
ecliptic pole to predictions based on ground-based observations and a stellar at- 
mosphere model. The brightness of point sources is calibrated to an accuracy of 
3%, relative to models for A star stellar atmospheres, for observations performed 
and analyzed in the same manner as the calibration stars. This includes correc- 
tions for location of the star in the array and the location if the centroid within 
the peak pixel. Long-term stability of the IRAC photometry was measured by 
monitoring the brightness of A dwarfs and K giants (near the north ecliptic pole) 
observed several times per month; the photometry is stable to 1.5% (rms) over 
a year. Intermediate-time-scale stability of the IRAC photometry was measured 
by monitoring at least one secondary calibrator (near the ecliptic plane) every 
12 hr while IRAC is in nominal operations; the intermediate-term photometry is 
stable with a 1% dispersion (rms). One of the secondary calibrators was found 
to have significantly time-variable (5%) mid-infrared emission, with period (7.4 
days) matching the optical light curve; it is possibly a Cepheid variable. 



1. Introduction 

The Infrared Array Camera (IRAC) was built at NASA Goddard Space Flight Center 
under the direction of the Smithsonian Astrophysical Observatory (Fazio et al. 2004). It 
is the mid-infrared camera on the Spitzer Space Telescope (Werner ct al. 2004), with four 
arrays observing at 3.6, 4.5, 5.8, and 8 /xm. The absolute calibration of the camera was 
performed in flight by comparing observed to predicted brightness for a set of stars that 
was selected and characterized before launch (Cohen 2003; Cohen et al. 2003). This paper 
presents the in-flight calibration results, including the observing strategy, the predictions and 
measurements, and an assessment of the calibration accuracy and stability of the instrument 
and the pipeline-processed data provided by the Spitzer Science Center to observers. 

2. Calibrator Selection and Observing Strategy 

The primary cahbrators (Table 1) are a set of A main sequence and K giant stars. 
They were chosen from a larger hst of candidates on the basis of having good spectral and 
photometric data and no known evidence of variability (Megeath et al. 2003). The set of 
plausible calibrators was observed during the in-orbit checkout to select which ones would 
become the primary calibrators for the nominal mission. Table 2 lists the coordinates of the 
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candidate primary calibrators that were observed during the in-orbit checkout, together with 
the IRAC channels in which useful data were obtained. Table 1 lists the primary calibrators 
that have been observed during the nominal mission. 

Primary calibrators are observed with all 4 IRAC arrays, every instrument campaign. 
The primary calibrators are located in the Spitzer continuous viewing zone, a region ap- 
proximately 10° in radius centered on the north ecliptic pole. The south ecliptic pole is also 
continuously visible, but it saves observing time to have all calibrators close together on the 
sky. After the first six months of the nominal science mission, the radius of the continuous 
viewing zone was decreased to 7?5, making two of the primary calibrators visible for only 8 
months each year. 

Secondary calibrators are observed approximately every 12 hr while IRAC is on; they 
are intended to monitor the gain stability during each observing campaign. Because routine 
Spitzer operations require pointing the high-gain antenna at the Earth every ~ 12 hr to 
downlink data, we chose a network of secondary calibrators that can be rapidly observed 
after each downlink. The Spitzer orbit is close to the echptic plane, and the telescope points 
opposite the high-gain antenna. Thus when the antenna points at Earth, the telescope points 
somewhere near the ecliptic plane at a longitude that changes with a period of about one 
year. Secondary calibrators spaced along the ecliptic can therefore be observed efficiently 
before or after telemetry downlinks, with little slewing overhead. The network of secondary 
calibrators has two stars every 20° of ecliptic longitude. One or two secondary calibrators 
are chosen for each campaign. 

All calibration star observations are performed by placing the star in 5 positions: the 
center of the array and the center of the four quadrants of the array. A single frame is taken 
at each position. Exposure times are set to place the peak brightness approximately 1/3 of 
saturation, so that the signal-to-noise is high and the nonlinearity corrections are accurate. 
The IRAC arrays saturate at ~ 10^ electrons or 40000 DN (whichever is lower). 

3. Photometry 

As input, we use the "basic calibrated data" (BCD) generated by the IRAC science 
pipeline (version SIO ^). In summary, the pipeline removes the electronic bias, subtracts 
a dark sky image generated from observations of relatively empty sky near the ecliptic 



^See the IRAC Data Handbook and pipeline history log and Pipeline Description Document at 
http://ssc.spitzer.caltech.edu/irac/data.html for details. 



Table 1. IRAC Primary absolute calibrators 



# 


Star'' 


Other name 


2MASSJ'' 


Type 


Av 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


1 


NPMlp64.0581 


HD 180609 


19124720+6410373 


AOV 


0.387 


9.117 


9.123 


9.104 


9.107 


9.122 


2 


HD 165459 




18023073+5837381 


AlV 


0.285 


6.584 


6.593 


6.575 


6.579 


6.591 


3 


NPMlp60.0581 


BD+60 1753 


17245227+6025508 


AlV 





9.645 


9.666 


9.653 


9.659 


9.669 


4 


1812095 




18120956+6329423 


A3V 


0.378 


11.596 


11.275 


11.254 


11.255 


11.265 


5 


KF08T3 




17551622+6610116 


K0.5III 





11.090 


10.944 


10.986 


10.985 


10.943 


6 


KF09T1 


GSC 04212-01074 


17592304+6602561 


KOIII 





8.114 


8.046 


8.087 


8.086 


8.047 


7 


KF06T1 




17575849+6652293 


K1.5III 





10.987 


10.762 


10.906 


10.833 


10.793 


8 


KF06T2 




17583798+6646522 


K1.5III 


0.189 


11.273 


11.065 


11.193 


11.133 


11.095 


9 


NPMlp66.0578 


GSC 04229-01455 


19253220+6647381 


KlIII 





8.330 


8.249 


8.347 


8.310 


8.335 


10 


NPMlp67.0536 


SAO 17718 


17585466+6747368 


K2III 


0.161 


6.409 


6.306 


6.428 


6.389 


6.348 


11 


NPMlp68.0422 


BD-l-68 1022 


18471980+6816448 


K2III 


0.462 


6.816 


6.712 


6.822 


6.780 


6.742 



''The star names used in this project are either abbreviated RA-only 2MASS designations, the designation NPMl from the 
Lick Northern Proper Motion program (Klemola et al. 1987) with "p" replacing +, the designation KF from a north ecliptic 
pole optical/near-infrared survey (Kiimmel & Wagner 2000). These are not intended to replace the formal designations. 

''Designation of the corresponding 2MASS catalog entry, with '2MASSJ' prefix removed, comprising the .12000 right ascension 
and declination sans punctuation (2 decimal places of seconds of right ascension and one decimal place of arcseconds of 
declination) . 
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Table 2: IRAC calibrator candidates 



Star" type RA Dec (J2000) Channels 



1732526 


A3V 


17 


32 


52, 


.6 


+71 


04 


43 


.1 


123 


1739431 


A3V 


17 


39 


43, 


.2 


+61 


02 


55, 


.8 


1234 


1740346 


A5V 


17 


40 


34, 


.7 


+65 


27 


15, 


.0 


1234 


1743045 


A5V 


17 


43 


04, 


.5 


+66 


55 


01, 


.7 


1234 


1757132 


A3V 


17 


57 


13 


.2 


+67 


03 


40, 


.0 


1234 


1802271 


A3V 


18 


02 


27, 


.2 


+60 


43 


35, 


.6 


1234 


1805292 


AlV 


18 


05 


29, 


.3 


+64 


27 


52, 


.1 


1234 


1808347 


A3V 


18 


08 


34, 


.8 


+69 


27 


28, 


.7 


1234 


1812095 


A2V 


18 


12 


09, 


.6 


+63 


29 


42, 


.3 


1234 


1812523 


AlV 


18 


12 


52, 


.4 


+60 


02 


32, 


.0 


12 4 


KF01T3 


Kim 


18 


03 


45, 


.5 


+66 


56 


03 


.0 


1234 


KF01T5 


Kim 


18 


04 


03, 


.8 


+66 


55 


43, 


.0 


1234 


KF02T1 


KOlll 


18 


02 


01, 


.7 


+66 


37 


39, 


.0 


1234 


KF02T3 


KOIII 


18 


03 


21, 


.8 


+66 


34 


08, 


.0 


1234 


KF03T1 


KOIII 


17 


57 


43, 


.9 


+66 


26 


55, 


.0 


1234 


KF03T2 


K1.5III 


17 


57 


51, 


.4 


+66 


31 


03, 


.0 


1234 


KF03T4 


Kim 


17 


59 


03, 


.9 


+66 


30 


59, 


.0 


1234 


KF05T1 


Kim 


17 


59 


12, 


.1 


+66 


41 


36 


.0 


1234 


KF06T1 


K1.5111 


17 


57 


58, 


.5 


+66 


52 


29, 


.0 


1234 


KF06T2 


K1.5111 


17 


58 


37, 


.9 


+66 


46 


52, 


.0 


1234 


KF06T3 


Kim 


17 


58 


50, 


.2 


+66 


49 


40, 


.0 


1234 


KF07T2 


Ki.5m 


18 


01 


40, 


.7 


+65 


59 


44, 


.0 


1234 


KF07T3 


K1.5III 


18 


02 


14, 


.8 


+66 


03 


00, 


.0 


1234 


KF08T3 


K0.5III 


17 


55 


16, 


.2 


+66 


10 


11 


.0 


1234 


KF09T1 


KOIII 


17 


59 


23, 


.0 


+66 


02 


56, 


.0 


1234 


NPMlp61.0569 


K0.5III 


17 


23 


25, 


.9 


+61 


12 


40, 


.7 


1234 


NPMlp61.0570 


Kim 


17 


26 


11, 


.6 


+61 


00 


48, 


.1 


1234 


NPMlp61.0582 


Kim 


17 


36 


55, 


.6 


+61 


40 


58, 


.1 


1234 


NPMlp64.0590 


KOlll 


19 


25 


14, 


.9 


+65 


01 


58, 


.9 


1234 


NPMlp66.0507 


K2111 


17 


31 


22, 


.1 


+66 


46 


35, 


.3 


1234 


NPMlp66.0513 


K2III 


17 


43 


53, 


.7 


+66 


31 


49, 


.8 


1234 


NPMlp66.0514 


K0.5III 


17 


45 


45, 


.4 


+66 


48 


45, 


.4 


1234 


NPMlp67.0536 


K2III 


17 


58 


54, 


.5 


+67 


47 


37, 


.4 


1234 


NPMlp68.0412 


K2III 


18 


28 


24, 


.8 


+68 


58 


51, 


.5 


1234 



" see note a of Table 1. 
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pole, flat-fields using a super-flat generated from the flrst year's calibration observations of 
relatively blank flelds near the ecliptic plane, and linearizes using laboratory measurements of 
each pixel's response to a calibration lamp in frames of varying length. For each BCD image 
in each standard star observation, aperture photometry was used to determine the source 
flux. The target is located as the image maximum after spatial filtering (low-pass, median, 
7 pixel width) to reduce noise and cosmic rays; this filtered image was used only for source 
identification. (Some of the candidate calibrators observed during in-orbit checkout have 
other stars nearby; they were located using the coordinates of the star and the astrometric 
calibration of the images.) The target location is then refined using a centroiding algorithm 
(cntrd in IDLphot). We place the image into electron units for proper error estimation; from 
the pipeline-processed images (which arc in MJy sr~^), we multiply by the gain (GAIN in the 
header; elec/DN) and exposure time (EXPTIME in the header) and divide by the calibration 
factor (FLUXCONV in the header; MJy sr~^ DN~^ s) that had been used in the pipeline; the 
scaled image at pixel i,j is then lij in electrons. To get the absolute brightness of each pixel, 
we also added the brightness of the zodiacal light at the time of the sky-dark observation 
(approximately 0.036, 0.18, 1.6, 4.4 MJy sr"-*^ in channels 1, 2, 3, and 4, respectively); this 
is only important for error propagation as it is a constant over the image. The background 
is calculated by a robust average in an annulus spanning 12 to 20 pixel radius centered on 
the target 

I sky — ^^sk^^^hj/Nsky (l) 

The on-source fiux is calculated by summing over a 10-pixel radius centered on the target 
aper in IDLphot). 

■P'elec ^on-^ij ^sky (2) 

The electron production rate, A^e = -^eiec/ 1 is proportional to the stellar fiux. 

An array-location-dependent photometric correction must be applied to IRAC photom- 
etry, to account for the variation in pixel solid angle (due to distortion) and the variation of 
the spectral response (due to the tilted filters and wide field of view) over the array (Hora 
et al. 2005). The photometric corrections are defined to be unity in the center of the array, 
so their role in the data analysis here is to remove a systematic error that would make the 
calibration star stars, which were observed at 5 widely-separated positions on the array (see 
Table 4 for the exact locations), have dispersions significantly higher than the measurement 
uncertainties. Table 4 shows the percentage corrections that were applied; specifically, if the 
table entries are fij then the fluxes in channel i at position j were divided by 1 + fij/ 100. 

A photometric correction must also be applied to IRAC photometry to account for the 
variation in the flux of a source as a function of location within a pixel. Sources placed 
directly in the center of a pixel give higher count rates than those placed near a pixel edge. 
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Table 3. IRAC Secondary calibrators for each campaign 



W 


Cm in 
ampaign 


Star 


RA (J2000) 


Dec 






61 




\A 51 


15 81 




18 01 


1 


IRAC003500 


HD 218528 


23 


08 


49.4 


-08 


48 


27 


KOIII 


6.311±0, 


,025 


6.367+0.045 


6.361+0, 


,024 


6.309+0.021 






SA115 554 


23 


41 


31.0 


+ 01 


26 


26 


K1.5III 


9.071±0 


,021 


9.213+0.025 


9.154+0, 


,018 


9.099+0.019 


2 


IRAC003600 


HD 4182 


00 


44 


25.4 


+09 


57 


09 


A2V 


8.990±0, 


,021 


8.981+0.021 


8.974+0, 


,021 


8.980+0.021 






SA116 554 


23 


41 


31.0 


+01 


26 


26 


K1.5III 


9.071±0, 


.021 


9.213+0.025 


9.154+0, 


.018 


9.099+0.019 


3 


IRAC003700 


BP 20 417 


02 


32 


14.4 


+20 


49 


09 


K1.5III 


7.777±0, 


.022 


8.058+0.025 


7.998+0, 


.018 


7.944+0.018 


4 


IRAC003800 


BP 14 649 


03 


18 


20.3 


+ 15 


25 


41 


K2III 


7.104±0. 


.021 


7.234+0.027 


7.192+0, 


.018 


7.138+0.019 


6 


IRAC003900 


HD 37725 


05 


41 


54.4 


+29 


17 


51 


A3V 


7.876±0, 


.023 


7.870+0.023 


7.864+0, 


.023 


7.869+0.023 






HD 244937 


05 


34 


54.2 


+30 


03 


28 


AlV 


8.512±0, 


,023 


8.488+0.023 


8.472+0, 


,023 


8.487+0.023 


6 


IRAC004000 


HD 55728 


07 


14 


47.6 


+ 16 


35 


50 


A2V 


9.465±0, 


,019 


9.464+0.019 


9.461+0, 


,020 


9.463+0.020 






HD 55677 


07 


14 


31.3 


+ 13 


51 


37 


A4V 


9.189±0, 


.020 


9.187+0.02 


9.184+0, 


.020 


9.187+0.020 


7 


IRAC004100 


AKS95 408 


08 


48 


52.6 


+ 19 


03 


42 


KOIII 


9.187±0, 


.020 


8.611+0.044 


8.606+0, 


.023 


8.463+0.019 






HD 77823 


09 


05 


15.8 


+ 14 


14 


43 


K1.5III 


6.183±0, 


.021 


6.325+0.026 


6.266+0, 


.018 


6.211+0.018 


8 


IRAC004200 


HD 89562 


10 


20 


08.9 


+03 


11 


56 


KOIII 


6.949±0, 


.025 


7.004+0.046 


6.998+0, 


.024 


6.946+0.021 






HD 90593 


10 


27 


44.2 


+04 


28 


20 


KOIII 


7.049+0, 


.026 


7.105+0.046 


7.099+0, 


.025 


7.047+0.023 


9 


IRAC004800 


HD 99253 


11 


25 


14.9 


-05 


52 


08 


K1.5III 


6.239±0, 


,022 


6.381+0.026 


6.322+0, 


,020 


6.267+0.020 


10 


IRAC004900 


HD 109441 


12 


34 


47.7 


-11 


46 


47 


KOIII 


7.123+0, 


,024 


7.179+0.044 


7.173+0, 


,024 


7.121+0.021 






HD 113745 


13 


05 


57.0 


-08 


41 


33 


AOV 


8.835±0, 


.020 


8.832+0.02 


8.829+0, 


.020 


8.834+0.020 


11 


IRAC005000 


HD 131769 


14 


55 


51.2 


-12 


12 


28 


K1.5III 


6.662±0, 


.022 


6.682+0.026 


6.616+0, 


.019 


6.667+0.020 






HD 135248 


15 


14 


21.4 


-08 


50 


39 


KlIII 


6.596±0, 


,025 


6.704+0.032 


6.663+0, 


,021 


6.611+0.022 


12 


IRAC005100 


HD 137429 


15 


26 


29.7 


-14 


28 


35 


KOIII 


6.181+0 


027 


6.227+0.046 


6.217+0 


,026 


6.169+0.023 






HD 145152 


16 


09 


51.2 


-12 


19 


09 


K2III 


6.570±0 


(124 


6.705+0.029 


6.666 + 


,022 


6.609 + 0.022 


13 


IRAC()()52()() 


HD 155334 


17 


12 


01.2 


-25 


40 


54 


AlV 


7.716±0 


,022 


7.714+0.022 


7.710+0 


(123 


7.714+0.023 






HD 156896 


17 


20 


55.5 


-20 


21 


35 


AOV 


9.364+0, 


,021 


9.355+0.021 


9.349 + 


,021 


9.357+0.021 


14 


IRAC005300 


HD 184837 


19 


37 


16.4 


-24 


40 


11 


AOV 


9.595±0, 


,021 


9.587+0.021 


9.581+0, 


,021 


9.588+0.021 


IB 


IRAC005400 


HD 195061 


20 


29 


33.2 


-18 


13 


12 


AOV 


9.330±0, 


.022 


9.318+0.022 


9.310+0, 


.022 


9.319+0.022 



Table 4: Array-location-dependent Photometric Corrections (%) 



index column 


row 


1 


Channel 
2 3 


4 


128 


128 


0.00 


0.00 


0.00 


0.00 


1 64 


192 


0.52 


0.53 


3.86 


-1.31 


2 192 


192 


0.46 


1.75 


-0.85 


1.50 


3 192 


64 


1.34 


1.56 


-1.63 


2.66 


4 64 


64 


1.34 


0.68 


3.39 


-0.39 
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The effect may be due to a nonuniform quantum efficiency across eadi pixel. We refer to the 
'pixel phase' as the distance of the centroid of a star from the center of the pixel containing 
the centroid. The pixcl-phase-corrections are normalized such that randomly-placed sources 
will have equal chances of requiring an increase or decrease in the flux. For randomly placed 
sources this location is 1/ ^J2t\ pixels from the center. The pixel phase dependence is only 
clearly detected in channel 1, and an approximate fit is 

Jvhase = 1 + 0.0535(l/v^ - V) (3) 

where p is the distance (in pixels) from the source centroid to the center of the pixel containing 
that centroid. Observed fluxes in channel 1 were divided by fphase- The channel 1 calibration 
stars tended on average to fall closer to the centers of their pixels than a random distribution, 
and the median correction (fphase) — 1-0%. 

The photometric uncertainty estimate includes three terms. First, the uncertainty in 
the subtracted sky background is the rms of the pixel values in the sky annulus, agky, divided 
by the square-root of the number of pixels in the sky annulus and multiplied by the number 
of on-target aperture pixels: 

^sky ^sky-^on/ 's/ -^sky i (4) 

where '^sky denotes the root-mean-square (rms) intensity in the sky annulus. Second, noise 
due to sky variations within the on-source aperture is 



sky, cm '^sky \/ -^on ■ (5) 

And third, the Poisson noise on the total background-subtracted counts within the on-source 
aperture is 

^poisson \/ ^elec- (6) 

The uncertainty in the flux is then the quadrature sum of the three error terms: 



ap — \ -I- Csfcu.on + '^loisson- i^) 



i,on ' poisson 

We have not explicitly included random variations in the instrument or its calibration; in- 
stead we attempt to measure those variations through the dispersion in photometry on a 
large number of images. All three terms in equation 7 contribute significantly to the pho- 
tometric uncertainty. For example, for NPMlp64.0581, the proportion (ysky'-'^sky,on'-<^poission 
is 0.64:0.66:0.39 for channel 1, 0.71:0.54:0.45 for channel 2, 0.69:0.58:0.43 for channel 3, and 
0.80:0.33,0.50 for channel 4. The differences in proportion are due to background brightness 
and dispersion, source flux, and exposure time in each channel. The Ugky can be reduced 
somewhat by increasing the size of the sky annulus. Both asky,on and asky can be reduced by 
decreasing the size of the on-source aperture, but then we would lose some flux and require 
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larger aperture corrections. The choice of on and sky apertures used here is optimized for 
cahbration purposes, where we attempt to gather as much of the star's flux as possible. For 
source detection, a smaller on-source aperture is recommended. All of the calibration stars 
are bright; for example, for NPMlp64.0581, the fiux-to-uncertainty ratio, Fjap is 158, 106, 
123, and 96 for channels 1, 2, 3, and 4, respectively, in a single BCD image. 

The flux and photometric uncertainty were calculated for each image of each calibration 
star; 15,341 fluxes were measured for IRAC campaigns 1-15. 

4. Photometric stability 

To monitor the stability of the camera on various timescales, the primary and secondary 
cahbrators were all reduced uniformly and plotted together. Figure 1 shows the results for 
the second IRAC nominal operations campaign. 

4.1. Stability on weekly timescales 

The relative repeatability of IRAC was measured by taking the root-mean-square varia- 
tions of each secondary calibration star photometry during a campaign. The distribution of 
these variances among all of the campaigns has a probability distribution function. The 
median of these rms values is 1.7%, 2.2%, 1.5%, and 1.6% in channels 1, 2, 3, and 4, re- 
spectively. The variance of the rms is consistent with the median photometric measurement 
error, which is 0.5%, 0.9%, 0.7%, and 0.9% in channels 1, 2, 3, and 4, respectively. Some 
stars have higher dispersions, but only one star at one wavelength (compared to a total of 
20 stars at 4 wavelengths) has an rms dispersion larger than 5%. This is a measure of the 
repeatability of photometry within a campaign, which can be summarized as being stable to 
2%. 



4.2. Stability on long timescales 

The camera is evidently stable on all timescales that have been checked. Table 5 sum- 
marizes the statistics characterizing the long-term stability for 8 primary calibrators over 11 
months. To assess the stability, we use the statistics for individual cahbrators in Table 5. 
Normalizing the count rates for each primary calibrator by their median, the long-term dis- 
persion in the photometry becomes 1.7%, 0.9%, 0.9%, and 0.5% in channels 1, 2, 3, and 4, 
respectively. 
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To determine whether the photometry has significant long-term variations (due to source 
or instrument), we inspected the long-term trends for each calibrator. Figure 6 shows the 
trend of the count rates over the first 11 months of the mission, and the last two columns in 
Table 5 show least-squares-fitted slopes to the temporal variations. An instrumental trend 
will be the same in the for each primary calibrator, while stellar variability will appear only in 
the individual star. Measurement errors generate dispersion (both random and systematic) 
that should not correlate between stars. Only one (KF09T1) of the 32 primary trend plots 
(8 calibrators, 4 channels) shows a slope that is formally significant. ^ The other primary 
cahbrators are all stable at a level of < 1.5%/yr. 



For each primary calibrator, we predicted the count rate in each IRAC channel. The 
procedure was explained in detail by Cohen et al. (2003). In brief, the spectral type of each 
star was determined from ground-based visible spectroscopy (Megeath et al. 2003). The 
intrinsic spectral energy distribution for that type was determined from models for the A 
stars (Kurucz 1993), or from empirical spectral templates derived from observations for the 
K giants (Cohen et al. 1996). For each star, an absolute normalization and the extinction Ay 
were derived by fitting the intrinsic SED to optical (UBV, Hipparcos Hp, Tycho-2 BTyr), 
near infrared (2MASS), and mid-infrared (IRAS, MSX) photometry. Table 1 shows the 
2MASS J, H, and Ks magnitudes (Skrutskie 1999) and the derived extinction for each 
primary calibrator. 

The reddened stellar SED template for each primary calibrator gives its fiux density at 
all wavelengths on a scale consistent with the network of standards constructed by Cohen et 
al. (1999) and absolutely validated by Price et al. (2004). The fiuxes in the IRAC wavebands 
were calculated by integrating the normalized spectral template for each star over the IRAC 
spectral response. The detailed derivation of the calibration factors is as follows. The number 
of electrons per unit time collected from a source with spectrum using a telescope with 
area A is 



5. 



Absolute calibration 



5.1. Flux prediction for primary calibrators 




(8) 



^This example shows how a simple monitoring program, with observations of 5 frames each, every 10-20 
days, can limit variability at the 2% level. Such a monitoring program takes approximately 2 hr of observing 
time per year, including slews to the target and observing at all 4 wavelengths. 
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Table 5: Long term (1 year) stability of IRAC primary calibration stars 



Sl.a.r 


(.'liaiiiic'l 


Akxliaii" 




\\ I'liicaii' 


WTbig'' 






NPMlp64.0581 


1 


0, 


.9672 


0.0163 





.9709 


0.0024 


-0.01246 


0.006787 


NPMlp64.0581 


2 


1, 


.0105 


0.0099 


1 


.0071 


0.0029 


0.00316 


0.004725 


NPMlp64.0581 


3 


1, 


.0080 


0.0082 





.9970 


0.0012 


0.00651 


0.003850 


NPMlp64.0581 


4 


1, 


.0017 


0.0125 


1, 


.0081 


0.0024 


0.00142 


0.005372 


HD165459 


1 





.9892 


0.0192 





.9815 


0.0004 


-0.01192 


0.008299 


HD165459 


2 


1, 


.0269 


0.0067 


1 


.0253 


0.0020 


-0.00212 


0.002973 


HD165459 


3 


1, 


.0412 


0.0083 


1 


.0422 


0.0008 


0.00366 


0.003180 


HD165459 


4 


1, 


.0390 


0.0053 


1 


.0373 


0.0012 


0.00127 


0.002308 


NPMlp60.0581 


1 


0, 


.9513 


0.0124 


0, 


.9497 


0.0018 


-0.00802 


0.007658 


NPMlp60.0581 


2 


0, 


.9931 


0.0183 


1 


.0205 


0.0008 


-0.01825 


0.010200 


NPMlp60.0581 


3 


1, 


.0018 


0.0828 





.9933 


0.0020 


-0.01570 


0.010688 


NPMlp60.0581 


4 


1, 


.0295 


0.0178 


1 


.0208 


0.0020 


0.00867 


0.011478 


KF09T1 


1 


1, 


.0268 


0.0082 


1 


.0292 


0.0010 


0.00806 


0.005187 


KF09T1 


2 


1, 


.0510 


0.0073 


1 


.0521 


0.0017 


0.00285 


0.004707 


KF09T1 


3 


1, 


.0195 


0.0069 


1 


.0178 


0.0011 


0.01244 


0.003559 


KF09T1 


4 


1, 


.0424 


0.0079 


1, 


.0470 


0.0011 


0.00365 


0.005214 


KF06T1 


1 


1, 


.0621 


0.0241 


1 


.0736 


0.0019 


-0.02734 


0.011175 


KF06T1 


2 


1, 


.0081 


0.0083 


1 


.0086 


0.0018 


-0.00551 


0.003902 


KF06T1 


3 


1, 


.0174 


0.0643 


1 


.0072 


0.0055 


0.04192 


0.024166 


KF06T1 


4 


1, 


.0648 


0.0323 


1 


.0745 


0.0048 


0.02743 


0.014229 


NPMlp66.0578 


1 


1, 


.0004 


0.0143 





.9971 


0.0004 


-0.00989 


0.006577 


NPMlp66.0578 


2 


0, 


.9934 


0.0096 





.9993 


0.0019 


-0.00448 


0.004144 


NPMlp66.0578 


3 





.9869 


0.0081 


0, 


.9853 


0.0013 


0.00651 


0.003476 


NPMlp66.0578 


4 





.9986 


0.0043 


1 


.0010 


0.0008 


-0.00018 


0.002175 


NPMlp67.0536 


1 


1, 


.0351 


0.0176 


1 


.0295 


0.0014 


0.00070 


0.007963 


NPMlp67.0536 


2 


1, 


.0507 


0.0097 


1 


.0558 


0.0025 


0.00085 


0.004416 


NPMlp67.0536 


3 


1, 


.0383 


0.0092 


1 


.0374 


0.0014 


0.00062 


0.004163 


NPMlp67.0536 


4 


1, 


.0470 


0.0044 


1 


.0463 


0.0010 


0.00253 


0.001915 


NPMlp68.0422 


1 


1, 


.0349 


0.0170 


1 


.0336 


0.0014 


-0.01374 


0.007078 


NPMlp68.0422 


2 


1, 


.0697 


0.0122 


1 


.0667 


0.0027 


-0.00912 


0.004984 


NPMlp68.0422 


3 


1, 


.0497 


0.0077 


1 


.0496 


0.0017 


0.00767 


0.003115 


NPMlp68.0422 


4 


1, 


.0605 


0.0057 


1 


.0597 


0.0010 


0.00180 


0.002322 



" Median of calibration factors measured in each of the first 11 campaigns, normalized to 
nominal values 0.1104, 0.1390, 0.6024, 0.2083 (MJy/sr)/(DN/s) for channels 1, 2, 3, and 
4, respectively. Root-mean-square of normalized calibration factors. Weighted mean of 
normalized calibration factors; with the weight for each point being the root-mean-square dis- 
persion of measurements during a campaign Formal statistical uncertainty in the weighted 
mean of the normalized calibration factors. Least-squares-fitted slope (and its statistical 
uncertainty) for the normalized calibration factors, in units of fractional variation per year. 
Apparently statistically significant slopes are in italics. 
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where R is the system spectral response (in units of electrons per photon, at frequency u, 
determined from pre-launch measurements). The calibration factor is the ratio of the flux 
density at the nominal wavelength, Aq = c/z/q, if the source had the nominal i/F,^=constant 
spectrum, to the observed electron production rate: 

C'^^. (9) 



where the color correction 



^ ^ j{F,/F,,:){v/vo)-'Rdv 



J {y/yo) ^ Rdv 



For an imager, like IRAC, we calibrate in units of surface brightness, and the flux measure- 
ments are made using aperture photometry in a finite aperture that does not necessarily 
include all of the flux. We will continue to use Fy as the true source spectrum, and define 
the calibration factor so that after both aperture and color corrections to the observations 
are performed, the best measurement of the source flux is obtained. To convert from raw 
IRAC units of data numbers (DN) into surface brightness per pixel, the calibration factor is: 

C C (11) 

where G is the gain (electrons/DN) and is the pixel solid angle (pixels are square with 
1.22" sides), and fap is the aperture correction factor (taken to be unity; see §5.5); the units 
of the calibration factor C are (MJy/sr)/(DN/s). The * superscript indicates the quantity 
refers to a calibration star, rather than a generic spectrum. It is perhaps worth emphasizing 
that the only properties of the stellar spectra that enter into the absolute calibration are 
the integrated in-band fluxes for each IRAC channel, and the cahbration factors are directly 
proportional to the in-band flux. 

Table 6 shows the predicted fluxes of the IRAC primary calibrators. The fiux predic- 
tions, F*^K*, are in IRAC 'quoted' units, such as would be obtained from science products 
from the SSC data analysis pipelines. The monochromatic flux at the nominal wavelength 
can be determined (for the calibrators or for any pipeline-processed IRAC data) by dividing 
the 'quoted' flux by the color correction (eq. 10). The nominal wavelengths for the IRAC 
bands were chosen to minimize the color corrections. One can show (by making a Taylor 
series expansion of the source spectrum) that the minimal dependence of the color correc- 
tion on the slope of the source spectrum is obtained for the nominal wavelength deflned as 
follows: 
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The nominal wavelengths calculated in this way, using the current spectral response curves 
(Fazio et al. 2004), are 3.550, 4.493, 5.731, and 7.872 fim in channels 1, 2, 3, and 4, respec- 
tively. Figure 7 shows the model spectra for an A dwarf and K giant calibrator, together 
with the spectral responses. We calculated the calibrator color-corrections K* separately for 
each calibrator spectrum. For reference, the color corrections are (in channels 1, 2, 3, and 
4, respectively) for an AlV star: K% =1.019, 1.001, 1.026, and 1.042; for an K1.5111 star: 
=1.021, 1.070, 1.001, and 1.093; for Rayleigh- Jeans spectra: Krj =1.011, 1.012, 1.016, 
and 1.034. 

The IRAC spectral response calibration convention, such that no color corrections are 
needed for sources with z/Fi,=constant spectra, is essentially the same convention used by 
IRAS (Beichman et al. 1988), COBE/DIRBE (Hauser et al. 1998), and ISO (Siebenmorgen 
et al. 1999), except that the monochromatic flux densities are at the nominal wavelength 
(eq. 12) for each channel rather than a round-number wavelength. 

5.2. Compctrison between different standard stars 

To search for systematic trends in the comparison between observed and predicted fluxes, 
we used two samples: the primary standard stars selected for the nominal mission (Table 1; 
11 stars) and the candidate cahbrators observed during the in-orbit checkout (Table 2; 34 
stars). The primary standards defined the IRAC calibration, so we consider the candidate 
calibrators as a relatively independent comparison sample to search for trends. Figure 8 
shows the histograms of the predicted/observed fluxes for the candidate calibrators, for each 
IRAC channel. The 'observed' fluxes have been scaled using the calibration factors averaged 
over the primary standards in the method described below in §5.3. The distributions are 
centered near unity, demonstrating good agreement between that the candidate calibrator 
and primary standards. 

We found one systematic trend in the predicted/observed fluxes: a displacement between 
the calibration derived from A dwarfs and K giants, in IRAC channels 1 and 2. Figure 8 
shows the histogram of the calibration factors separately for these two types of primary 
cahbrators. The difference in average calibration inferred from the two types of cahbrators 
is twice as large than the widths (3.4%) of the distributions for each calibrator type. This 
leads to an unnecessary systematic error in the calibration, with a spectral dependence that 
changes channels 1 and 2 with respect to the other IRAC channels. To bolster the statistics, 
we include the primary standards together with the calibrators observed only in IOC; the 
total sample comprises 13 A stars and 29 K stars. The difference between the A and K star 
distributions of predicted/observed fluxes is 7.3 ±2.3% in channel 1, 6.5 ±2.4% in channel 2, 
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3.6 ± 2.5% in channel 3, and 2.1 ± 2.8% in channel 4, with the predicted/observed higher for 
K stars than for A stars in all four channels. Our original calibration strategy was to observe 
both A and K calibrators, then take the average to guard against model uncertainties. After 
seeing such a clear systematic difference, with a color dependence, between the two types of 
calibrators, and furthermore seeing that there is no temporal drift in the calibration (so that 
we can use all calibration star data all the way back to the in-orbit checkout), we decided to 
adopt an A-star-only absolute calibration convention. At present, we do not know why the 
K giants are systematically offset from the A stars, nor do we conclusively know which of the 
two types of calibration are correct. With the relatively small number of stars observed and 
analyzed here, and only two types of calibrators analyzed, it is not possible to conclusively 
say whether the template spectra are inaccurate representations of the stars, or the particular 
sample of stars is anomalous. 

The distribution of predicted/observed flux (Figure 8) is non-gaussian, because it is a 
combination of the measurement uncertainties (which do have a gaussian distribution) and 
the systematic uncertainties (which do not). In all four IRAC channels (but particularly 
channel 1), two peaks in the distribution are evident due to the separate distributions of 
A and K stars. Further, there are too few calibrators for the central limit theorem to lead 
to a clean gaussian distribution for each type. The dispersions of the distributions are 
~ 2.5 — 3.5%, which is comparable to the uncertainty in the predicted fluxes based only on 
pre-launch ground and space measurements with other telescopes. 



5.3. Best calibration factors and their uncertainties 

The recommended calibration factors for IRAC were calculated as a weighted average of 
calibration factors for the four. A- type primary calibrators. The results are given in Table 7^. 
Weights were inversely proportional to the uncertainties in the absolute calibration of each 
star (Cohen et al. 2003). 

The total uncertainty in the calibration factor is not the statistical uncertainty in the 
weighted mean, since aabs = 1-5% of the uncertainty in each of the predicted fluxes arises 
from the absolute cahbration of Vega and Sirius, which scales all of the predictions. To 
assess the uncertainty in the cahbration, we separate three items in Table 7. The absolute 
calibration uncertainty for each star is (t ground- The dispersion in calibration factors de- 
rived from an individual star is (Jrms-, attributed largely to measurement errors. Finally, the 
campaign-averaged calibration factor for each star has each calibrator has some dispersion 



^The calibration factors in Table 7 have been applied to IRAC data with pipeline version Sll. 
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from campaign to campaign of arepeat, which hmits long-term drifts (which are apparently 
not present for IRAC). The combined uncertainty for n calibrators determined from 

= ^J (Tibs + (<o«nd - ^lbs)h' + ^rms/n- (13) 

We used this equation, with n = 4, to derive the uncertainty (2.0% in all four IRAC bands) 
in the absolute calibration factors listed in Table 7. 

Including the larger sample of calibrators observed only during IOC, the calibration 
differs from the values in Table 7 by -0.9, -1.1, -1.7, and -2.3 %, respectively, in channels 1, 
2, 3, and 4. This difference is marginally significant and we incorporate it into the overall 
calibration uncertainty. The IRAC calibration will be updated in the future to include these 
results and a wider sample of spectral types from the Legacy science programs. 



5.4. IRAC magnitude system 

We define the IRAC magnitude system such that an observer measures the flux density, 
pguot^ of a source from the calibrated images out of the IRAC pipeline, performs image- 
based corrections (array-location-depcndcnt photometric correction, pixel-phase correction, 
aperture correction), and then uses the zero-magnitude flux densities, Fjero, to calculate the 
magnitude \i] = 2.5 logio(i^iio/i='r"*) where i = 3.6,4.5,5.8,8 are the four IRAC channels. 
In this system, there is no need to know the spectral shape of the source, as the magnitude 
is a measure of the in-band flux relative to that of Vega. The zero-magnitude fluxes were 
determined by integrating the Kurucz model spectrum of a Lyr over the passbands using 
the equations above; speciflcally. 

The resulting zero-magnitude flux densities are 280.9 ± 4.1, 179.7 ± 2.6, 115.0 ± 1.7, and 
64.13 ± 0.94 Jy in the [3.6], [4.5], [5.8], and [8] //m channels, respectively. Due to the choice 
of A dwarfs as the absolute cahbrators, this magnitude convention should yield results on 
the same scale as used in optical and near-infrared astronomy. 



5.5. Aperture corrections and extended emission 

All the discussion so far has been based on aperture photometry with a specifled beam 
size and sky annulus (§3), and the calibration factors in Table 7 are deflned for these pa- 
rameters. The parameters were optimized for photometry of isolated, bright point sources 
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Table 6. Predicted brightness of Primary absolute calibrators 



Predicted flux, F^Jf* (mjy)°- 



# 


Star 


3.6 /im 


4.5 //m 


5.8 jum 


8.0 (um 


Campaign log 


1 


NPMlp64.0581 


63.03 ± 1.49 


41.03 ± 0.95 


26.17 ±0.62 


14.40 ± 0.34 


yyyyyynnnny 


2 


HD165459 


647.38 ± 17.02 


421.19 ± 10.85 


268.71 ± 7.06 


148.09 ± 3.90 


yyyyyyyyyyy 


3 


NPMlp60.0581 


38.20 ± 0.89 


24.74 ± 0.56 


15.75 ±0.37 


8.70 ±0.20 


nnnnyyyyyyy 


4 


1812095 


8.68 ± 0.20 


5.66 ± 0.12 


3.62 ±0.08 


2.00 ± 0.04 


yyyyyyyyyyy 


5 


KF08T3 


11.78 ± 0.42 


7.25 ± 0.45 


4.64 ±0.13 


2.69 ± 0.08 


yyyyyyyyyyy 


6 


KF09T1 


169.94 ± 5.86 


104.69 ± 6.33 


67.02 ± 1.74 


38.76 ± 1.04 


ynnnyyyyyyy 


7 


KF06T1 


13.93 ± 0.46 


7.80 ± 0.50 


5.34 ±0.13 


3.09 ±0.08 


yyyyyynnnny 


8 


KF06T2 


10.53 ± 0.35 


5.99 ± 0.38 


4.05 ±0.10 


2.34 ± 0.06 


nnnnyyyyyyy 


9 


NPMlp66.0578 


140.94 ± 6.11 


82.39 ± 7.34 


54.54 ± 1.60 


29.72 ± 1.03 


yyyyyynnnny 


10 


NPMlp67.0536 


843.49 ± 35.09 


482.35 ± 39.32 


319.84 ±9.81 


185.25 ± 7.39 


yyyyyyyyyyy 


11 


NPMlp68.0422 


580.54 ± 24.89 


335.64 ± 27.73 


223.18 ±7.01 


128.96 ± 5.22 


yyyyyyyyyyy 



^Predicted flux density in 'quoted' IRAC convention. The predicted flux density at the nominal wavelength is F*^ , and the 
color correction for each star's predicted spectrum spectrum is K*. 



Table 7: Absolute calibration factors for IRAC" 



Channel 


Calibration factor" 


(Aground 


^rms 


'^repeat 


1 


0.1088 ±0.0022 


2.3% 


2.0% 


0.9% 


2 


0.1388 ±0.0027 


2.3% 


1.9% 


0.9% 


3 


0.5952 ±0.0121 


2.3% 


2.1% 


0.9% 


4 


0.2021 ± 0.0041 


2.3% 


2.1% 


0.8% 



units: (MJy/sr)/(DN/s) 
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and will not in general be suitable for other applications. In particular, some of the source 
flux lands in the sky annulus and is subtracted, while additional flux lands outside the sky 
annulus and is ignored. Measurements using different beam sizes or sky annuli will need to 
account for these effects. Therefore, 

1. Point sources extracted with the same aperture as we used for calibration will get the 
correct flux (no further aperture correction needed); 

2. Point sources extracted with other on-source radii or sky annuli can use aperture 
correction tables, as long as they are normalized such that there is no correction for 
apertures identical to those used in our absolute calibration; 

3. Point sources extracted using point spread function fitting should verify the normal- 
ization of the point spread function would give unity flux if the on-source and sky 
calibration apertures are applied (i.e. not unity flux for the integral over the entire 
point spread function); and 

4. extended emission surface brightness will be incorrectly calibrated and will require 
scaling by the "inflnite" aperture correction. 

The reason we have not attempted an aperture correction is that the measurements are 
still under way and the empirically-derived values are not well understood. Part of the effect 
is simple diffraction and could be estimated using, for example, an Airy function determined 
by the primary mirror size and the flltcr central wavelength; this effect explains channels 1 
and 2. However, in channels 3 and 4 substantially more of the source flux is scattered out of 
the calibration aperture than can be explained by diffraction theory. This light is thought 
to be scattered within the detector material. Its distribution, dependence on location of the 
source in the array, and the ultimate fate of the lost source flux are under investigation. 
The empirical determination of the aperture corrections is being described by Marengo et 
al. (2005); at present we recommend extended emission (including the basic calibrated data 
and the image mosaics from the pipeline) be multiplied by the effective aperture correction 
factors of 0.944, 0.937, 0.772, and 0.737 for channels 1, 2, 3, and 4, respectively. 

Three further practical matters have already been mentioned but are worth summarizing 
here as well. 

1. The IRAC science pipeline generates images in units of surface brightness, but because 
of distortion, the pixels do not subtend constant solid angle. 

2. The IRAC spectral response varies over the field of view, and therefore the color cor- 
rections are field dependent. 



-18- 



3. The electron rate in the 3.6 /xm channel depends slightly on pixel phase. The present 
calibration factors is correct on average, as is appropriate for sources observed multiple 
times at multiple dither positions. For the most precise photometry, however, pixel 
phase should be taken into account. 

Corrections (1) and (2) are available from the Spitzer Science Center and are described by 
Hora et al. (2005); they were taken into account in the present paper by using the array- 
location-dependent photometric corrections in Table 4. Correction (3) was applied using the 
simple equations in §3. Observers should apply these corrections in a manner consistent with 
that applied in this paper, in order for the absolute calibration to apply. 

6. Serendipitous results 

The secondary calibrator SA115_554 has a periodic light curve. It was observed every 
12 hr in two campaigns, spanning a total of 26 days. The light-curve phases well to a period 
of 7.4 days, with an amplitude of 5%. Figure 9 shows the hght curve, phased to at modified 
Julian date 52996.5. The light curve is identical at all 4 IRAC wavelengths, indicating that 
the source of variability is not a mid-infrared-specific phenomenon but rather something 
with colors comparable to the star, which is classified as a K1.5III. This star was found to 
be variable as part of the AS AS by Pojmanski (2000), who found I-band variability with a 
period of 7.451 days, identical to the IRAC period. The amplitude of the I-band variation 
was 0.18 mag, and its hght-curve is distinctly different from that measured by IRAC. If 
the variability were due to an orbiting companion, then the companion would probably be 
late-type and small, perhaps an M dwarf. But the observed variability likely arises from 
pulsations of the star itself. The optical light curve is more characteristic of a Cepheid 
variable, and the spectral type obtained by Drilling & Landolt (1979) was G7, consistent 
with a Cepheid, though the spectral type derived from our own work was a K1.5III. 

The primary calibrator HD 165459 was found to have a significant (~ 40%) excess at 
24 //m, in the MIPS calibration program. Presumably, this excess is due to a disk around 
the star such as in other 'Vega'-type A stars. Is there any evidence for this disk at IRAC 
wavelengths? Calibrating the measured fiuxes with the new calibration factors, we find the 
star's fiux relative to the photosphcric model is 0.982±0.020, 0.987±0.024, 0.976±0.022, and 
0.993±0.022, in channels 1, 2, 3, and 4, respectively. Since these are absolute measurements, 
we have used arms from Table 7 as the calibration uncertainty (and combined with the 
uncertainty in this individual star's measurements). In fact, we can remove the absolute 
uncertainty in the size of the star by normalizing all fluxes at 3.6 /im. Then the ratio 
of observed flux to photosphcric, assuming there is no excess at 3.6 /xm, is 1.006 ± 0.009, 
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0.994± 0.009, and 1.011 ±0.008 in channels 2, 3, and 4, respectively. It appears that there is 
no mid-infrared excess for this star in the 3.6-8 fim range, at a 95% confidence upper limit 
of 4% (with no assumptions), or at 4.5, 5.8, or 8 /im, at a 95% confidence limit of 1.2, 1.8, 
or 1.6% (assuming no excess at 3.6 /im). 

7. Conclusions 

The Infrared Array Camera on the Spitzer Space Telescope has a stable gain on all 
measured timescales, making it possible to measure variability at the 2% level for carefully- 
reduced data. The absolute calibration using stellar photospheric models scaled to ground- 
based photometry at optical through near-infrared wavelengths is accurate to 1.8%, 1.9%, 
2.0%, and 2.1% in channels 1 (3.6 /um), 2 (4.5 /xm), 3 (5.8 /im), and 4 (8 /um), respectively. 
To measure fluxes at this high level of accuracy requires several photometric corrections: 
array position dependence (due to changing spectral response and pixel solid angle over 
the camera of view), pixel phase dependence (due to nonuniform quantum efficiency over a 
pixel), color correction (due to the different system response integrated over the passband 
for sources of different color) , and aperture correction (due to the fractions of light included 
within the measurement aperture and lost in the background aperture). 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through an award issued by 
JPL/Caltech. This publication makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Massachusetts and the Infrared Processing 
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Fig. 1. — Summary of calibration star results during the second IRAC campaign of the 
routine science mission. Each point is a measurement of the absolute calibration factor, 
normalized by the new cahbration factors derived in this paper [0.1104, 0.1390, 0.6024, 
0.2083 (MJy/sr)/(DN/s)]. The clusters of points and the beginning of each plot are the 
primary calibrators, which always have the same symbol type in each plot. Filled symbols 
are K giants and open symbols are A dwarfs. The primary calibrators (and their symbols) 
are NPMlp64.0581 (filled square), HD 165459 (filled triangle pointing up), NPMlp60.0581 
(filled diamond), 1812095 (filled triangle pointing down), KF08T3 (open square), KF09T1 
(open triangle pointing up), KF06T1 (open diamond), KF06T2 (open triangle pointing 
down), NPMlp66.0578 (open circle), NPMlp67.0536 (open triangle pointing right), and 
NPMlp68.0422 (open triangle pointing left). Two error bars are shown for each point, indi- 
cating the rms dispersion and the uncertainty in the mean. Spread throughout each campaign 
are the secondary cahbrators, observed every 12 hr. One or more secondary cahbrators are 
observed each campaign, with source names indicated in a legend for each campaign. One 
of the two secondary calibrators (SA115-554) is apparently a variable star and has been 
removed from our calibrator list. 
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Fig. 2. — Same as Figure 1, but for campaign IRAC004100. 
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Fig. 3. — Same as Figure 1, but for campaign IRAC004900. 
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Fig. 4.— 



Same as Figure 1, but for campaign IRAC005200. 
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Fig. 5. — Same as Figure 1, but for campaign IRAC005400. This campaign length is typical 
of the second year of routine operations. 
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Fig. 6. — Long-term monitoring of the calibration factor derived from individual primary 
calibrators. Each panel is for a different IRAC channel: 1 (top left), 2 (top right), 3 (bottom 
left), and 4 (bottom right). Within each panel, the calibration factors for each primary 
calibrator have been normalized to their median, and they have been offset from the previous 
curve by 0.1. 
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Fig. 7. — Predicted spectra of two primary calibrators, an A dwarf and a K giant, together 
with the relative spectra responses of the four IRAC channels. Both the spectra (flux density 
units; oc Jy) and relative spectral response (in electrons per photon) were scaled arbitrarily 
to unity to fit on this plot. 
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Fig. 8. — Distributions of predicted/observed fluxes for the 34 calibrator candidates in each 
IRAC channel. Each box represents one star. The histogram of the total sample is the 
upper envelope. The histogram for A stars only has diagonal hatching one way (/), and the 
histogram for K stars only has diagonal hatching another way (\). The two distributions 
overlap near unity, but the displacement is evident, with the K stars having higher predicted 
than observed flux. 
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Fig. 9. — Light curve of the secondary cahbrator SA115_554 in the four IRAC channels. The 
sohd curve (same exact curve in each of the four panels) is a fourth-order polynomial fit to 
all of the wavelengths, indicating the light curve is color- invariant. 
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